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Abstract. Neutral particles can be guided and focussed using electric field 
gradients that focus in one transverse direction and defocus in the other, 
alternating between the two directions. Such a guide is suitable for transporting 
particles that are attracted to strong electric fields, which cannot be guided using 
static fields. Particles are only transmitted if their initial positions and transverse 
speeds lie within the guide's phase space acceptance. Nonlinear forces are always 
present in the guide and can severely reduce this acceptance. We consider the 
effects of the two most important nonlinear forces, a term in the force that is 
cubic in the off-axis displacement, and a nonlinear term which couples together 
the two transverse motions. We use approximate analytical techniques, along 
with numerical methods, to calculate the influence of these nonlinear forces on 
the particle trajectories and on the phase space acceptance. The cubic term alters 
the focussing and defocussing powers, leading either to an increase or a decrease 
of the acceptance depending on its sign. We find an approximate analytical result 
for the phase space acceptance including this cubic term. Using a perturbation 
method we show how the coupling term leads to slow changes in the amplitudes 
of the transverse oscillations. This term reduces the acceptance when it reduces 
the focussing power, but has little influence when it increases that power. It is 
not possible to eliminate both nonlinear terms, but one can be made small at the 
expense of the other. We show how to choose the guide parameters so that the 
acceptance is optimized. 
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1. Introduction 

A neutral particle in an electric field gradient feels a force. The field polarizes the 
particle and, if it has a gradient, pulls on the induced dipole moment. In its ground 
state, a neutral particle is always attracted to regions of strong electric field. In 
other quantum states, the attraction may either be to strong fields, or to weak fields. 
Here, we are concerned with the guiding and focussing of strong- field-seeking particles. 
Static electric fields cannot do the job and so a dynamic focussing scheme needs to be 
used [nil]- An alternating gradient (AG) focusser is a series offenses, each focussing in 
one transverse direction and defocussing in the other. The focus and defocus directions 
alternate. Particles have stable trajectories through the focusser because they are close 
to the axis when they travel through the defocussing lenses, and further away when 
passing through the focussing ones. 

The purpose of an alternating gradient guide is to transmit a large number 
of particles over arbitrarily long distances with as little loss as possible. Such a 
guide might employ an electrode structure that is uniform along the beamline, with 
the alternation achieved by switching the applied voltages - alternation in time. 
Alternatively, the voltages can be fixed with the alternating gradients set up by the 
electrode structure itself - alternation in space. Early experiments with alternating 
gradient structures used two-rod and four-rod geometries to focus polar molecules, 
e.g. [2111]. These same ideas were later applied to cold atoms, demonstrating that 
an atomic beam or fountain could be focussed this way [SJ [Hj. Of considerable 
interest at present is the application of AG focussing to the Stark deceleration of 
strong-field seeking molecules. Prototype alternating gradient decelerators have been 
demonstrated experimentally [II [3 [9], and their focussing properties investigated 
in detail both theoretically and experimentally [10]. For heavy molecules to be 
decelerated to rest, these decelerators will need to employ very many stages of AG 
focussing, typically 100 lenses or more. Recently, a long AG guide with a double bend 
has been used to guide slow molecules from an effusive source [TT]. As buffer gas 
cooling technology advances, it is likely that such guides will be more widely used as 
efficient tools for extracting molecules from the cold buffer gas [121 IB] • An AG guide 
has also been used to separate conformers of C6H7NO by exploiting the selectivity 
of the guide to the ratio of dipole moment to mass [14] . The theoretical analysis we 
present here is suitable for describing all of these linear guiding structures. It can also 
be applied directly to other alternating gradient structures such as a storage ring for 
strong-field seeking molecules ^15^ or linear ac traps of the kind demonstrated for both 
molecules [TS] and atoms [17]. The physics we present also carries over directly to ac 
traps with cylindrical symmetry because the same terms appear in the equations of 
motion [18] . 

The anharmonic forces present in AG guides and traps for neutral particles tend 
to be detrimental to the phase space acceptance, often reducing this by a factor of 
10 or more. Similarly, these nonlinearities reduce the quality of the images in an AG 
focussing system. This problem has long been recognized. In reference [1] for example, 
molecules were guided in a four rod geometry with the rod positions carefully chosen 
to optimize the linearity of the forces, resulting in a large improvement in transmission 
over previous two rod geometries. The optimization of electrode geometry with 
respect to minimizing nonlinear forces has been discussed in [19j and |10j . Numerical 
simulations have been used to predict the performance of several types of AG guides 
[TT] , decelerators [lOj and traps [TS] [16] all confirming the hugely deleterious effects of 
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the nonlinearities. 

In this paper we use a mixture of analytical and numerical techniques to calculate 
the effect of the nonlinear forces on the particle trajectories and the phase space 
acceptance of an AG guide. We start by reviewing the linear theory since this is the 
foundation for the analysis of nonlinearities. It is well known that the motion consists 
of an oscillation driven at the frequency of alternation, called the 'micromotion', 
superposed on a slower oscillation of larger amplitude called the 'macromotion'. 
A good approximate description of the macromotion is obtained using an effective 
potential. The anharmonic terms modify this effective potential and so alter the 
phase space acceptance. The leading anharmonic terms in this potential are a quartic 
term along with a term that couples together motion in the two transverse directions. 
The quartic term in the potential leads to terms in the force that are cubic in the 
transverse coordinates. We treat this quartic term analytically, deriving an expression 
for the phase-space acceptance in the presence of this term. To understand the effect 
of the transverse coupling term we apply a perturbation method which reveals the 
nature of the coupled trajectories and shows that this term tends to decrease the 
phase space acceptance even when it increases the depth of the effective potential. We 
use numerical techniques throughout to support the theory, test the approximations 
and obtain results when the perturbation approach ceases to be valid. We apply our 
analysis to some example AG guides for neutrals. 

2. Equations of motion 

Our starting point is equation (8) of reference [1 Oj . which gives a suitable form for 
the electrostatic potential inside an AG guide, $(a;, y), as a function of the transverse 
coordinates x and y: 



This equation is derived by writing down a multipole expansion of the potential, and 
then eliminating terms that are unsuitable for a guide. The scale factors, ro and $0 
characterize the size of the electrode structure and the applied voltages, while ai, 03 
and as are dimensionless constants representing the sizes of the dipole, hexapole and 
decapole terms in the expansion. 

A neutral particle inside the guide will be polarized by the electric field to a degree 
that depends on the particle polarizability and on the magnitude of the field. The 
interaction of this induced polarization with the electric field changes the potential 
energy of the particle; this is the Stark shift, W . If the electric field is spatially 
inhomogeneous, so too is W , and there is an associated force acting on the particle, 
F = — VVF. For our present purpose, we do not need to know the detailed dependence 
of the Stark shift on \E\, the magnitude of the electric field. We need only know that, 
for sufficiently small electric fields, the Stark shift is quadratic in \E\, W = — ia|i?p, 
while for sufficiently large electric fields it becomes linear, W — —ii\E\. The constants 
of proportionality in the two cases are known as the polarizabilty, a, and the dipole 
moment, ^. Whether the electric field is considered small or large depends on the 
polarizability of the particle. For example, the Stark shift of a CaF molecule is in the 
linear regime at a field of lOOkV/cm, whereas the Stark shift of a Li atom is quadratic 
in the same field. 




(1) 
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We now impose the additional constraint that as <C 03 <C ai, meaning that the 
electric field in the guide is dominated by a constant term, and that the decapole term 
is much smaller than the hexapole term. We then obtain the following approximate 
expression for the n'th power of the electric field magnitude, valid throughout the 
region where x,y < tq: 



\E{X, = E^{1 + asniX^ - Y^) + (a^n(3 - n) - 6a5n)X'^Y^ 

+ {aln{n - l)/2 + a5n){X^ + y*) + • • •)• (2) 

In this equation, we have set the redundant scale parameter ai to unity and have 
introduced scaled transverse coordinates X — x/ro, Y — y/r^. Eq \s the electric field 
at the centre of the guide. The neglected terms are of order 0305, a\ and higher. Since 
n can be any number, we see from ^ that the potential in which the particles move 
always has the same functional dependence on the transverse coordinates, irrespective 
of whether the Stark shift is linear, quadratic, or intermediate. Although we have 
treated the case of a long guide, similar terms appear in the treatment of cylindrically 
symmetric traps [18], and so our analysis could also be readily applied to that case. 

In writing the equations of motion, we may either use the axial coordinate, z, 
or the time, t, as the independent variable. The former is most appropriate for long 
guiding structures, while the latter tends to be used for trapping geometries. We will 
use z throughout, writing the equations of motion as 



X"{Z) = 77^ y^X + xY'X + C^' 

Y"{Z) = + xX^Y + CY^^ , (3) 

where the primes denote differentiation with respect to the scaled axial coordinate 
Z = z/r^^ and 77, x s-^d C ^-re dimensionless parameters that depend on the character 
of the Stark shift, but not on the spatial coordinates. For a linear Stark shift, 
X = 2|a3| - 605/103!, C = Sas/laal, and 

For a quadratic Stark shift, x— kal ~ 6a5/|a3|, ^ = laal + 2a5/|a3|, and 

We have introduced the quantity 7?., which characterizes the Stark shift on the axis of 
the guide in terms of the forward kinetic energy, and will be useful later in the paper. 
We have written the equations of motion for a lens that focusses in the y-direction 
when 03 is positive. For focussing in the a;-direction, simply interchange X and Y . 

We note from ([3]) that the dominant term in the force is linear in the transverse 
coordinate, focussing in one plane and defocussing in the orthogonal plane. In addition 
to the desirable linear term, there is a term cubic in the transverse displacement, and 
also a nonlinear transverse coupling term. By a suitable choice of the parameters 05 
and a3, it is possible to eliminate one or other of the nonlinear terms, but simultaneous 
elimination of both nonlinear terms is impossible without also setting 77 to zero. For 
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example, when the Stark shift is hnear, the cubic term in the force can (in principle) 
be eliminated by choosing 05 = 0; alternatively, the coupling term can be removed by 
choosing 05 = a§/3; the elimination of both terms is clearly not possible. Thus, the 
lenses in an alternating gradient guide for neutral particles are necessarily aberrant. 
More generally, the values of C, and x ^-re linearly related via x + 3C = 2n|a3|. This 
relationship shows that it is also impossible for both x C to be negative, implying, 
as we will show later, that no matter how the electrodes are arranged at least one of 
the terms is detrimental to the transmission of the guide. 

In ([3]), we are considering the case where 05 ^ 03 ^ 1, and so the coefficient 
of the linear term is much larger than the coefficients of the nonlinear terms. It is 
therefore tempting to suppose that the effect of the nonlinearity on the trajectories is 
small. However, in going from a focussing lens to a defocussing lens, the linear term 
changes sign (as it must) but the nonlinear terms do not. The alternating gradient 
guide can support a stable beam because the envelope of the beam is always smaller in 
the defocussing lenses than in the focussing lenses. For this reason, the net effect of the 
linear forces is to confine the molecules, but this trajectory-averaged confining force 
is considerably weaker than the linear forces of the individual lenses. By contrast, the 
nonlinear forces are identical in the two lenses and simply add up. It follows that, 
even though the nonlinear terms in the individual lenses are small, their role is crucial 
in understanding the dynamics of the guided particles. 



3. Motion in the linear approximation 

Throughout this paper we consider a linear AG guide consisting of alternate converging 
and diverging lenses each of length L, separated by gaps of length S. To understand 
the effect of the nonlinear forces, we first analyze the motion in the absence of the 
nonlinearities. This idealized motion has been discussed in some detail in reference 
[To] and is based on the original work of Courant and Snyder who discussed the 
same problem in the context of charged particle focussing. We give a brief summary 
of this theory here, and extend the discussion of reference [10) . 



3.1. Trajectories 

Neglecting the nonlinear terms, we write the equation of motion in one dimension as 

X"{Z) + K{Z)X{Z) = Q. (6) 

Here, K{Z) — ri'^G{Z), where G(Z) is a function whose value is +1 in a focussing lens, 
— 1 in a defocussing lens, and in a drift space. We know the solution in each region 
of constant G. For example, the solution inside a focussing lens is given by 



X \ _ f cos{riZ) Tj-^sm{riZ) \ f X 

X' ) - I -rjsinivZ) cos{7jZ) ) [ X' 



Zo 



= F{r^,Z)(^A (7) 

V / Zq 

The solution inside a defocussing lens is identical to ^ , with the replacement of 77 by 
IT], and the transfer matrix denoted by D{r], Z). For a region of free space, let 77 ^ 
in ([7|) and call the matrix 0{Z). 
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Knowing the solutions in each of the separate regions, we can construct the fuU 
trajectories piecewise, but the Courant-Snyder formahsm offers a far better approach. 
This involves looking for a general solution of the form 



X{Z) ^ y/eP{Z) cos{iP{Z) + S) 



= Aiy/P{Z)cosij{Z)+A2^P{Z)simP{Z) (8) 

where /3 is a Z-dependent amplitude function that has the same periodicity as the 
AG array, "0 is a Z-dependent phase, and e, S, Ai and A2 are defined by the initial 
conditions. 

By substitution wc find that (H]) is a valid solution to ([6|) provided 

^P(Z) = f -^dZ (9) 
' Jo P{Z) 

and 

-\p'^ + \pP" + KP^ = l. (10) 

It appears that all we have achieved is the replacement of the differential equation for 
X{Z) by a more complicated differential equation for (3{Z). The advantage is that we 
will never have to solve this new differential equation. We already know the solution 
to (O, but piecewise, and our aim is to express this in a convenient form. From ([5]) 
we have 

X\Z) ~ ^^^{—acosip — sin?/>) + ^^(— asin^ + cos^') (11) 

where 

a = -1/3'. (12) 

Using ([5]) and (|lip we find the relationship between the coordinates X, X' at position 
Z and those at position Z + LccXh -^ceii being the periodicity of the AG array. We 
make use of the periodicity constraint on /3, (}{Z + Lceii) = fi{Z) (and hence also on 
a) and thus obtain: 



X \ I cos 11 + a sin fi (3 sin /i \ ( X 

X' / \ —J sin 11 cos /i — a sin /i J \ X' 



(13) 



where 7 = (1 + a )/ (i and 



li = V^iccii) = / -TTJ-dZ (14) 

Jz p{Z) 

is the phase advance per unit cell. Since the integral is taken over a full period, /i 
is independent of Z . The matrix appearing in (jl3p is known as the Courant-Snyder 
matrix, M.. 

The procedure to obtain j3{Z) is now straightforward. We illustrate this procedure 
for a simple guide with no gaps. Consider an alternating sequence of focussing and 
defocussing lenses, each of length L, with Z — defined to be at the entrance of a 
focussing lens. We first obtain an explicit expression for the transfer matrix between 
a point Z inside a focussing lens and the equivalent point one lattice unit further 




Figure 1. ID phase-space acceptance, in units of "qro, as a function of the 
parameter r}L. SoHd line: Uniform aperture throughout the guide. Dashed line: 
Aperture at defocussing lenses only. Dotted line: Approximate result for small 
micromotion. 



downstream. This matrix is M = F{j],L — Z).D{r], L).F{j], Z). We then calculate 
i Tr(M), which, from the Courant-Snydcr matrix is equal to cos/i. The explicit result 
is cos/x = cos{riL) cosh{r]L). Finally, we calculate (3 by equating the upper right hand 
element of M to that of Ai , obtaining 

cosh(r]L)sm(r]L) + cos(riL -2riZ)smh{riL) 
I3{Z) = , , (15) 

■qj 1 — cos^(r/L) cosh^(?7_L) 

for Q < Z < L. We can apply the same procedure to a point Z inside a defocussing 
lens (i.e. L < Z < 2L). The result for P{Z) is identical to ([15]) with the replacement 
of 77 by ir] and the replacement of Z hj Z — L. Since /3 has period 2L, we now know it 
everywhere, and since ipi^) can be calculated directly from (3{Z) we have the complete 
general solution to the equation of motion in the form of ([8|). This same procedure 
can equally well be applied to more complicated guiding structures. 



3.2. Phase-space acceptance 

Returning to the case of a general guide, particles enter with a range of transverse 
positions and angles and we wish to know which of them are transmitted. Consider a 
particle whose initial conditions are given by e and S, as in Using the first line of 
^ to form the quantity X'^ + {aX + PX'Y , and recalling that (3' = -2a, i}}' = 11(3 
and 1 + = 7/?, we find the invariant 

+ 2aXX' + (3X'^ = e. (16) 

Suppose we record the transverse position and angle of the particle at the 
longitudinal positions Z + NLceii, for a large number of integer values of N . Then 
a plot of these points on a phase-space diagram, whose axes are X and X', traces 
out the ellipse defined by (jl6|) . If we repeat this procedure for a different value of 
the shape of the ellipse evolves periodically with Z according to but its area 
remains the same, 7re. A distribution of particles having ti < e fills the ellipse. If 
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the guide has a constant transverse aperture with walls at X — ±1, then it is clear 
from ([5]) that all particles having < e will be transmitted indefinitely provided 
that \/ (3max^ < 1, where Pmax is the maximum value of f3{Z). If this condition is 
satisfied, the transverse extent of the beam never exceeds the transverse aperture of 
the guide. Particles that have > l/f3max may also be transmitted by the guide 
provided that the turning points of the macromotion (the cosine factor in ([5])) never 
coincide with the turning points of the micromotion (the amplitude function in ([S])). 
That condition may occur for a large number of particles if fi/ir is a rational number, 
/i being the phase advance per unit cell defined by (HH). However, these "resonances" 
in the phase-space acceptance become very narrowly centred on these special values of 
jj, once the guide is long compared to the wavelength of the macromotion, and we can 
usually ignore them. Then, the transverse phase-space acceptance is simply Tr/Pmax- 
As discussed in [10] , the maximum value of f3 always occurs in the centre of a focussing 
lens. Returning to our specific example of a gapless guide, this is explicitly seen in 
at Z = L/2. It follows that, in this case, the transverse phase space-space acceptance 
in one dimension is 



Since we are using Z as the independent variable, the acceptance has the units of 
XX', i.e. tq ■ radians. There are, of course, two transverse directions, and the 2D 
transverse acceptance is simply the square of the above formula. 

So far, we have supposed that the transverse aperture formed by the electrodes 
of the guide is uniform along its entire length. For many electrode geometries, this 
is a poor approximation. An interesting case is that of alternating two-rod lenses, 
discussed in |10j and used in all experimental work on alternating gradient deceleration 
reported to date. Here, the electrodes lie in the defocussing plane and there is no 
aperture at all in the focussing plane. The maximum beam envelope still occurs at 
the centre of the focussing lens, but there is no aperture at this point for the beam to 
crash into. Instead, the appropriate limiting value of /? to insert into the acceptance 
equation is its value at the entrance to the defocussing lens. As a result, the phase- 
space acceptance increases. For the case of a gapless guide with an aperture at the 
defocussing lenses only, the second term in the denominator of P7)) is replaced by 
cos{t]L) smh{r]L). The acceptance in one transverse direction of a gapless guide is 
plotted in figure \T\ for both kinds of transverse aperture. We see that when the 
aperture is defined only by the defocussing lenses (dashed line) the ID acceptance 
peaks at a slightly higher value of rjL where it is almost twice as large. Later in 
this paper, we will make a "small micromotion approximation" which amounts to 
neglecting the modulation of (3, replacing it by the constant value Lccii//^- lu this 
approximation, the phase space acceptance in one dimension becomes tt/x/LccH- For 
a gapless guide, this approximate acceptance is shown by the dotted line in figure [1] 
We see that the approximate result agrees well with the exact results when tjL <C 1, 
because the micromotion amplitude is indeed small in this regime. As r]L increases, 
the approximate result overestimates the acceptance. At r/L — 1 for example, the 
approximate acceptance is 31% larger than the exact value with uniform aperture, 
and 2% larger for an aperture at the defocussing lenses only. 

We turn now to some specific examples. Consider first a beam of ground state 
CaF molecules travelling at 350 m/s into a gapless guide having Eq — lOOkV/cm 




(17) 



cosh(?7L) sin(ry_L) + siiih(riL) 
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Figure 2. The fill factor as a function of the distance D (in units of tq) between 
the entrance of the guide and an upstream square aperture of side tq. The result 
shown is for the particular case of a gapless guide with r) = 0.07 and L = 20. 

and 03 = 0.2. In this case, 77 = 0.053, and with tq = 1 mm the acceptance reaches 
its maximum value of 39mmmrad when the lenses are 24 mm in length. Thus, the 
acceptance of the guide, which can transport the molecules over an arbitrarily long 
distance, is the same of that of a 1 mm^ square cross-section pipe of length 26 mm. As 
a second example, consider Li atoms launched with a speed of 10 m/s from a magneto- 
optical trap (MOT) into the same guide. Here, 77 = 0.31, the optimal lens length is 
4 mm and the transverse acceptance is 230mm.mrad, equivalent to 2.3mni.ni/s in 
position-velocity space. This acceptance is larger than the phase-space area occupied 
by the atoms in a typical Li MOT, indicating that all the atoms can be successfully 
transported. 

3.3. Fill factor 

The previous section showed how to calculate the transverse phase-space acceptance, 
but often the source is unable to fill this acceptance area completely. The final output 
of particles depends on the overlap between the acceptance and the phase-space area 
occupied by the particles at the guide's entrance, which we call the fill factor. As 
an illustration, consider a beam formed in a supersonic expansion, passing through a 
small aperture (a skimmer) and then into the guide. For the purposes of calculation 
it is convenient to proceed in the opposite direction. We calculate the phase-space 
acceptance at the entrance of the guide, project this back through the skimmer and 
onto the source to find out what fraction is filled by the source. In practice, the 
source is usually large enough and divergent enough to completely fill the part of the 
acceptance that fits through the skimmer, and we assume this to be the case in the 
calculation that follows. 

Let us write the phase-space ellipse accepted by the guide as 

AiX"^ + 2BiXX' + = 1. (18) 

Here, Ai = l3maxl{Q), Bi = /3,„aa;a(0) and Ci = l3maxfi{Q), where Z = is defined 
to be at the entrance of the guide. We project this ellipse through a distance —D 
and so obtain a new ellipse with coefficients A2, 82,02- The relationship between the 
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coefficients is 

( M\ / I Q \ ( Ai\ 

(19) 

This is easily derived using the drift matrix 0{Z) (see (O). 

At this point, we can calculate the area of that part of the phase-space ellipse 
that fits through the skimmer aperture. The aperture provides a constraint in the X 
direction, but not in the X' direction. Integrating the ellipse between X = ±R, we 
obtain the area [J, 

J [Bl - A2C2)X^ + C2 dX. (20) 

Repeating this for the other transverse direction, and multiplying together the two 
results yields the 2D transverse phase-space acceptance of the entire guide and 
aperture setup. This result, divided by the 2D acceptance of the guide alone, 
(T^/Pmax)'^, gives the fill factor. Figure [2] shows an example of how the fill factor 
decreases as D increases, stressing the importance of minimizing the distance between 
the guide and any aperture placed upstream. 



4. An effective potential approach to the nonlinear dynamics 

In order to solve the nonlinear problem, whose solution is X{Z), let us define new 
coordinates, w and s, in terms of the function /3(z) which appears in the solution to 
the linear problem, equation 



X{Z) 

' ^ - [ ^^^ = ^/'^' (21) 

The phase advance s increases by 1 each time the particle advances through one 
unit cell. The quantity w represents the macromotion of the particle. Applying the 
transformation to ^ we obtain 



X" + KX ^ 



+ (1/2 P (3" -a' + Kf3')f,'w 



^202 ' " . ^^r ,r - -0, (22) 

where the primes still denote differentiation with respect to Z. Note that we have 
made no assumptions about the form of X{Z) in this transformation. With the help 
of Uni) and (fT^ we see that 1/2 /3/3" — + Kp'^ ~ 1, and so our equation of motion 
in the transformed coordinates reduces to that of a harmonic oscillator of angular 
frequency /i, 

w + ^'^w = 0, (23) 

where the dots denote differentiation with respect to s. In the transformed coordinates 
the micromotion has disappeared and the macromotion appears as motion in a 
harmonic potential whose curvature is proportional to /i^. 



I This procedure is only accurate for a rectangular aperture, since we have assumed throughout that 
we can treat the two transverse directions independently, multiplying the results together at the end. 
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Re-introducing the nonlinear terms, we write ([3]) in the form 
X" + Tf{GX ~ G^C,X^ - G'^xY^X) = 0, 

Y" + rfi-GY - G^C,Y^ - G^xX^Y) - 0, (24) 

where we use the function G{Z) introduced earlier so as to handle an arbitrary 
sequence of converging lenses, diverging lenses and drift spaces. Applying the same 
coordinate transformation as above, this becomes 

wi + y.^{wi - G^CifPlwl - G^x-rfPlhwlwi) = 0, (25a) 

W2 + fi^{w2 - G'^CifPlwl - G^xv^Pif3iwlw2) = 0. (256) 

In pSap . wi, Pi and si are identical to the w, (3 and s introduced earlier, and the dot 
indicates differentiation with respect to si. The quantities W2, /?2 and S2 in (j256p are 
the corresponding quantities for the y motion, and differentiation is with respect to 
S2- Naturally, (ii and (52 are related: 132{Z + Lcou/2) = (ii{Z). 

So far, our analysis has been exact. In order to make progress, let us assume that 
the guide is operated well within the range of stability. In this regime the amplitude of 
the micromotion is small compared to the macromotion amplitude and the envelope 
function, (3{Z), is approximately constant. We set (3i — (32 = /3o- It follows from (fT4| 
and (j2ip that /3o = Lcc1i/m and Si — S2 — Z/Lcdi- In addition, the wavelength of the 
macromotion is large compared to the length of a unit cell in this regime, and so G^ 
oscillates rapidly relative to the macromotion wavelength. We will replace it by its 
mean value (G^). In this way, the same formalism will apply to guides with or without 
gaps. For many practical guides these approximations will turn out to be useful and 
not very restrictive. 

Making these replacements in (|25ap and (|25bl) , and then transforming back to the 
more familiar variables, X and Y, we obtain 

X" + ^^X - ebiX^ - eb2Y^X = 0, (26) 

where differentiation is with respect to s, interchange of X and Y yields the second 
of the two equations of motion, and 



eh=C{vLcen)^{G^}, (27) 
eb2=x{vLccii)^G^). (28) 

Patience will reward the reader's curiosity about the seemingly superfluous little 
epsilon. 

In the small micromotion approximation, we can think of the particles moving in 
an effective potential, Vcs. From ([26]) . and the equation obtained by interchange of X 
and Y, we find this effective potential to be given by 

l/eff/A/ = ^^^2(^2 ^ - \ebi{X^ + Y^) - \eb2X^Y\ (29) 

We also note that the transverse phase space acceptance obtained by applying the 
small micromotion approximation in the harmonic case is simply 7r/i. Converting 
from our present coordinates, (X, ^), back to (x, ^), this acceptance becomes 
TTfxrjro/lrjLccn)- This result was already discussed below (|17p in the context of figure 

m 
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5. The cubic term 

In this section we derive analytical expressions for the trajectories of the particles and 
for the phase space acceptance of the guide, with the cubic term included. We set 
62 = in (j26p . multiply by X' , and then integrate with respect to s, obtaining 

+ \^?X'^ - \ebiX* ^ h, (30) 

where ft, is a constant of the motion (the total energy) , and we identify the first term as 
the kinetic energy and the second and third terms as the potential energy. Rearranging 
for X' and integrating gives 

' ' (31) 



where p = [i} jlh and q = ebi/Ah. To complete the solution, we factorize the expression 
under the square root to the form (1 — a_|_X^)(l — a_X^) where a± = p/2± \/p'^/4: — q. 
Then, the integral takes the standard form for an elliptic integral of the first kind, F, 
and the solution i^ 

s{X) = -=L= F{^m-\^X), a^/a+). (32) 
^2a+ft 

Finally, this expression can be inverted to give X(s) in terms of a Jacobi elliptic 
function: 

1 



X{s) = sn(-\/2a+ft s, a_/a+). (33) 



This equation provides an analytical form for the trajectory of a particle through the 
guide, in the presence of the cubic nonlinearity, but averaged over the micromotion. 
Figure [3l^a) shows some examples of how the trajectories depend on the size of 
the cubic term, calculated using ([55]) . The figure shows that as ebi increases, the 
wavelength increases and the trajectories evolve from a sinusoidal shape towards a 
more square shape. These changes reflect the way the effective potential changes with 
ebi . Part (b) of the figure shows the effective potential along the X-axis for the same 
parameters used to obtain the trajectories plotted in (a). As the value of ebi increases, 
the steepness of the effective potential well decreases, and so the wavelength of the 
macromotion increases. In the extreme case where ebi — Z^^, the particle will stop at 
the maximum of the potential and remain there forever. 

From (l3T|) and the discussion that follows it, we find the wavelength of the 
macromotion to be 

4 dX 

V2ft7o ^1-a+XVl-a-^^ 

Xinax being the turning point of the motion where X' = 0, found to be X,„ax = 1/ 
The wavelength evaluates to 

A=^^A'(a_/a+), (35) 

§ Conventions for elliptic integrals are not standardized. We use the convention where F{4),m) = 
(1 — m sin^ (d))~^ d6. The inverse of this elliptic integral is the Jacobi amplitude am(u,m) 
such that if u = F{<p,m), then if> = am{u,m). The Jacobi elliptic function is defined by 
sn(u,m) = sin(am(u, m)). 
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Figure 3. (a) The macromotion for the case = 7r/4 and Xmax = 1 and for 
three different values of tb\. Solid line: tb\ = 0. Dotted line: eb\ = 0.25. Dashed 
line: ebi = 0.5. (b) Plots of the effective potential along the X axis, when b2 = 0, 
fi = Tt/A, and eb\ has the same three values as in (a). The equation for the 
effective potential is 1291 1. (c) Macromotion spatial frequency as a function of 
jfmax, when tb\ = 0.5. (d) Macromotion spatial frequency as a function of ebi, 
when Xmax =0.5. In (c) and (d) we have used the exact expression for v given 
by 1351 1. normalized to fo = /x/(27r). 



where K is the complete elUptic integral of the first kind. It is illuminating to write 
down the series expansion for A, or its inverse, the spatial frequency v, in powers of 
the small quantity tbi. 

A 27r V 8/^2 256^*4 

For positive ebi , the spatial frequency decreases with increasing ebi and with increasing 
oscillation amplitude. To lowest order, the change depends on the square of the 
oscillation amplitude. Particles that have a large amplitude will oscillate more slowly 
through the guide than those that remain close the guide's axis. Figure [3ljc) and (d) 
illustrate the exact spatial frequency versus X^ax and ebi. These plots shows the 
quadratic dependence on the oscillation amplitude and the linear dependence on ebi. 

We now calculate the phase-space acceptance of the guide in the presence of the 
cubic term. The two transverse directions are uncoupled at this point in our discussion, 
so it is convenient to calculate the acceptance in ID, squaring the result to get the 
total 2D acceptance. When ebi is positive, the effective potential in which the particles 
move has a turning point at X = Xc where ^^Xc — ^biX^ = 0, i.e. 
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Particles with enough energy to reach this turning point cannot be confined in the 
guide. A particle that has just enough energy to reach the top of the efi^ective potential 
well has an energy = 1/2^^X^ — = /i*/(4e6i). To determine the phase- 

space acceptance we need to distinguish between two cases. In the first case, the 
turning point is 'inside' the physical aperture and so the acceptance is limited by the 
turning point. In the second case, the turning point is 'outside' the physical aperture, 
or there is no turning point at all (because efoi < 0), and the acceptance is limited by 
the physical aperture. In the first case, the equation defining the area in phase space 
where particles stay inside the potential is 

+ _ igfe^x* = ///(4e6i). (38) 



The enclosed area is the phase space acceptance and is easily calculated by rearranging 
for X' and integrating under the curve from X = 0to X = X^ where the curve crosses 
the X-axis. The result is 

where we have used equation (|57|) . The phase-space acceptance scales as the cube of 
the phase-advance per unit cell, a measure of the trapping strength, and is inversely 
proportional to ebi, the coefficient of the nonlinearity. 

In the second case, where the acceptance is limited by the physical aperture at 
X — Xa , the equation defining the accepted phase-space area is 

+ i^^X^ - iebiX^ = i^^^l - \eb,X\. (40) 

By integrating the expression for X' with respect to X from X = to X = Xa, we 
find the area enclosed by this curve to be 



A = -^X^Vl - l/2{XA/X,y{E{K) - [1 - {Xa/X,)^]K{k)}, (41) 

where k — X\/{2X'^ — X'^), and K and E are the complete elliptic integrals of the 
first and second kinds. We note that, for Xa = 1, this equation reduces to A — tt^ 
in the limit where C ~* 0. This is the result we already found when applying the 
constant beta approximation in the harmonic limit. We find it useful to factor out 
the TT/i, writing our result as the product of the acceptance in the harmonic limit, Aq, 
and a factor, Fi^, that accounts for the cubic nonlinearity: 

A = AoF^, (42a) 

where 

F, = ^X^ (42&) 
when Xc < Xq, and 

F^ = Ax,V1-1/2(Xa/X,)2{£;(^) _ [1 - iXA/X,r]K{^)} (42c) 

OTT 

when Xc > Xq. By writing the result in this way, we can choose to use the exact result 
for Aq, derived using the procedure detailed in section [3.21 rather than 7r/i which we 
already know from figure [T] is not very accurate unless rjL is small. In this way, we 
eflFectively use the small micromotion approximation to handle only the nonlinear part 
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Figure 4. The effect of the cubic term, characterized by C,, on the ID acceptance 
of a gapless guide. Points and hnes are numerical and analytical results for (a) 
riL = 0.5, and (b) r?L = 1. 

of the problem. This approach, which we apply throughout, is justified by the high 
accuracy of the results it gives when compared to numerical simulations. 

Together, (j42&P and (|42c[) give the phase-space acceptance for all values of efci. 
The accuracy of this analytical approach depends on the validity of the approximation 
we have made, namely that [3 is approximately constant. To explore this, we can 
calculate the acceptance numerically for some specific cases, and then see how well 
the numerical and analytical results compare. To calculate the acceptance exactly, 
we set X = in (j24p and then solve this equation numerically for a large number 
of particles, taking care to ensure that the initial phase-space distribution is large 
enough to completely fill the acceptance area. Particles whose trajectories exceed the 
transverse boundaries of the guide (at X = ±1), are removed from the calculation 
so that we are left with the set of particles transmitted by the guide. We calculated 
for a guide consisting of 200 lenses. This was long enough to show that the phase- 
space acceptance is independent of the guide's length, provided this is longer than the 
wavelength of the macromotion. 

Figure [4] shows how the acceptance of a gapless guide depends on the value of (. 
It compares the exact acceptance, calculated numerically, with the analytical result, 
Fi^, given by (142 6p and (|42cp . Part (a) shows the results obtained with rjL = 0.5 
for which the phase advance per cell is /i = 0.14. The numerical results have been 
normalized to the exact acceptance obtained when C = 0, which in this case is 0A26r]rQ, 
while the analytical result is already normalized. When C is positive, the acceptance 
decreases rapidly with increasing as the turning point of the effective potential 
moves inwards. For negative values of C, the aberration increases the curvature of 
the effective potential and so increases the acceptance. The analytical and numerical 
results show excellent agreement for all values of C calculated. For positive ( the results 
agree to within the error bars on the numerical values, which are approximately the 
size of the points. For negative values of C, the analytical result is consistently 2% 
larger than the numerical results. Part (b) of the figure shows the results obtained 
with rjL = 1 for which fi — 0.58. Again, the numerical results have been normalized 
to the aberration-free acceptance, which is 0.701?7ro. Here, our small micromotion 
approximation begins to break down, and the true acceptance is less well represented 
by the analytical result. Nevertheless, the analytical result differs by 20% at most from 
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the true value. The analytical approximation is very useful since it can be evaluated 
instantly, unlike the numerical simulation. The full 2D acceptance is simply the square 
of the ID acceptance that we've calculated. 

We now examine the effect of the cubic term on some example guides. At present, 
all the experimental work on alternating gradient deceleration of polar molecules has 
used two-rod lenses. In these decelerators, two rods, 6 mm in diameter, with centres 
separated by 8 mm, are aligned with their axes parallel to the molecular beamline, 
and the experiments are performed in the linear Stark shift regime. In this case 
05/03 = 03 = 1/7 and so C, = 2/7. Usually, the length, S, of the drift spaces is 
approximately equal to the length of the lenses, and so we take r]S = rjL. We use 
our analytical formulae to calculate how the cubic term reduces the 2D phase space 
acceptance in this geometry, normalizing the results to the values obtained when C = 0. 
When TjL — 0.5, the normalized acceptance, Fq, is a miserable 3%. This grows to 18% 
when TjL = 0.75 and to 60% when rjL — 1. We see that the nonlinearity is hugely 
detrimental if the guide is operated at low values of rjL where the phase advance per 
cell, and hence the depth of the effective potential, is small. As rjL is increased towards 
the edge of stability, the nonlinearity becomes less important. This geometry is even 
worse for the electric guiding of cold atoms whose Stark shift is quadratic, because 
the quadratic Stark shift gives C = Sjasj instead of 2|a3|. 

In the four rod guide considered in [TU], rods of radius 2.3ro were placed with 
centres on the corners of a square of side 6.6ro. The rods on one side of the square 
were at positive high voltage, and those on the other side at negative high voltage, 
so the arrangement has two-pole symmetry. In this arrangement the cubic term is 
small when the Stark shift is linear, ( — 0.028, and the normalized acceptance, F(^, is 
20% when rjL = 0.5 and 79% when 77L = 1. The guide is not so suitable when the 
Stark shift is quadratic, for then ( = 0.17. The structure used in [11] to guide slow 
molecules from an effusive source was built from four rods of radius ro with centres on 
the corners of a square of side Srg. The rods on two diagonally opposite corners were 
at high voltage, ±V, with the other two grounded. Fitting the electrostatic potential 
obtained in this geometry to the multipole expansion, ([T]), we obtain 03 = 0.59 and 
05/03 = —0.055. So this geometry has ( — —0.11 for polar molecules with linear 
Stark shifts, meaning that the cubic term enhances the acceptance over the perfect 
linear case. However, the large value of 03 means that our series expansion of the 
electric field, equation ([2]), converges very slowly. The expansion does not provide an 
accurate representation of the field unless terms of higher order are included and so 
our conclusion about the effect of the cubic term is not so helpful in this case. We 
will see later that the lowest order transverse coupling term is severely detrimental for 
this geometry. 

6. Multiple scales analysis 

In the preceding section we found analytical solutions for the motion in the presence 
of the cubic nonlinearity. We can no longer do so when 62 in equation ([26]). 
Nevertheless, provided the nonlinear terms are small, we can elucidate the motion of 
the particles using a perturbation technique known as multiple scales analysis. We 
are particularly concerned with stable motion in the guide, i.e. with those particles 
that remain inside the guide indefinitely. We recognize that the coupling term will 
lead to an exchange of energy between the x and y oscillations, and so the oscillation 
amplitudes will be modulated. If e&2 is small enough, this modulation will occur on a 
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much longer scale than the wavelength of the macromotion. To handle this, we look 
for a solution in the form of a perturbation expansion using two separate length scales. 
We follow references [HI [551 H3] j and refer the reader to these books for more details 
on the multiple scales method. 

The independent variable in ([26)1 is s = Z/LccW- We introduce a second length 
scale, a = es, treat X and Y as functions of both s and a, and then expand X and 
F as a power series in e: X{s,(j) = Xo{s,a) + eXi{s,a) + . . .. In order that the 
initial conditions on X and Y be satisfied uniquely in this expansion, we choose to set 
^„(0) = r„(0) = X'JO) = y;^(0) = O for aU n O. The derivatives are 



dX dXa fdXo dXi 



ds ds \ da ds 
d'^X d^Xo f^d^Xo d^Xi 



ds"^ ds"^ \ dsda ds"^ 

Substituting into (j26p we obtain a differential equation whose solution should not 
depend on the arbitrary quantity e. We therefore have a separate equation for each 
power of e, these being 



ds 



2 



m'^o = 0, (43) 



^ + M^Xi = -2^ + h,X^, + h^Y^Xo, (44) 

and similarly for Yo, ^i, by simple interchange of the symbols X and Y . 
The solutions for Xq and Fq are 

Yo - ^e*(*'«+»2) + c.c, (45) 

where c.c stands for the complex conjugate and J^i , ?'2 1 j ^2 are real functions of tr, 
but independent of s. These solutions are now substituted into (j44p . giving, after a 
little algebra. 



5^ 

ds'^ 



fi'^Xi = Cle'^"' + cse^'"' + c.c, 



- I in'^''^ + nr. "^^^ + ^^V^ + 2\ , &2 2 .(2e2-9i) 



Consider for a moment the solution to the above differential equation with the 
right hand side replaced by zero (the corresponding homogeneous equation). Its 
solution has terms in e^'^" and e"*'^", terms which also appear in the right hand side 
of the inhomogeneous equation, l|46p . Terms in an inhomogeneous equation which are 
themselves solutions to the associated homogeneous equation lead to secular terms in 
the solution. These secular terms grow with s more rapidly than the corresponding 
solution of the homogeneous equation by at least a factor of s [23], i.e. as scos(/Lts) 
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in this case. If our solution is to represent stable motion, corresponding to particles 
trapped in the guide indefinitely, we must have to eliminate the secular terms by setting 
to zero the coefficient, Ci, of e^^" in pS)) . We are free to do this, since the r's and 0's 
are otherwise undetermined functions of a. Setting ci = and then separating out 
the real and imaginary parts, we obtain the following differential equations describing 
the slow variation of ri,r2,0i,O2: 

ri^ -^rlrl sm{2{e2~0i)) = O, (47a) 
da 

+ ^rirj sin(2(02 - ^i)) = 0, (47&) 

? + + ^r,'[2 + cos(2(02 - ^i))] = 0. (47rf) 

da Hfi 8/i 

By solving this set of differential equations we obtain the slow evolution of the 
amplitudes and phases, which we can then substitute into equations (05]) to get the full 
solution to lowest order in e. As we will see later, this procedure is remarkably accurate. 
A great deal of insight into the motion of the particles can be acquired from the above 
set of equations. By adding together (|47ap and (|476l) we see that + r| = Eq, a 
constant, proportional to the total energy at this order of approximation. Subtracting 
(|Tra|) from (|476p and P7c|) from (|47(i|) yields 



^-^i?oC(l-Osm(27)-0, (48) 

g + W-20(|i-|(2 + cos(2,))).0, (49) 

where ^ is the fractional energy associated with oscillations in the x direction, 
^ = rf/Eo, and j = 62 — Oi is the phase difference between the two transverse 
oscillations. The solutions of these equations tell us how the phase difference and 
the partitioning of the energy evolve with the axial coordinate. Even more instructive 
is to eliminate the independent variable by dividing ([49]) by (|48p. and then integrating 
to give 

{36i -b2[2 + cos(27)]K(? - 1) - C, (50) 

where C is a constant. This equation is the central result of this analysis. It tells 
us that there is an invariant quantity involving only ^ and 7. For a given particle, 
the value of C is fixed by the the initial conditions. The particle is then confined to 
travel on the contour in (^, 7) space defined by the above equation. This tells us how 
the amplitudes of the two transverse oscillations will change, allowing us to determine 
whether a particle will remain in the guide or crash into its boundaries. 
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Figure 5. Trajectories tlirougli an alternating gradient guide with nonlinear 
transverse coupling. The guide has riL = 1, rjS = 0, x = 0.03. The initial 
conditions of the particle were X(0) = 0.5, X'(0) = 0, y(0) = 0.25, y'(0) = 0.25. 
Solid line: Exact trajectory calculated numerically. Dashed line: Approximate 
trajectory given by I I45I I and the solutions to I l47all - M7d\ . The Y-trajectory has 
been offset (by 1 unit) for clarity. 



7. The nonlinear coupling 

7.1. Small coupling coefficient 

In this section, we use the results obtained above to understand the effect of the 
nonlinear coupHng on the motion of particles in the guide. To isolate the effect of 
the coupling, we set bi = 0. To demonstrate that the approximation developed 
above is a useful one, we begin by looking at an example trajectory. The solid 
line in Figure [5] shows the exact trajectory of a particle travelling through 120 unit 
cells of a gapless guide that has ryL = 1, C = and x — 0.03. The trajectory 
shows the usual superposition of micromotion and macromotion. The transverse 
coupling is responsible for the slow modulation of amplitude and relative phase 
evident in the figure. The dashed line in the figure shows the result of the multiple 
scales approximation. The parameters given above correspond to /i = 0.585 and 
£^2 — {'nLcc\\)^X = 0.12, and the dashed curve is obtained from ([IS]) using the solutions 
to (|47ap - ()47<ij) . We see that the approximate solution does a good job of describing 
the true motion of the particle, including the modulation of the amplitude and phase. 
In fact, in order that the approximate and exact trajectories could be distinguished 
in the figure, we had to choose a case where our approximations begin to break down 
- the micromotion amplitude is not so small relative to the macromotion, and £62 is 
not very small compared to /Lt^. In situations where the approximations hold more 
accurately, the two trajectories are virtually indistinguishable. 

Now that we have illustrated the usefulness of this approximate approach, we 
continue with our analysis, returning to the key result given in (|50p . When bi = 0, 
this equation reduces to 



representing a set of contours of a form that is independent of 62 (except as an over- 
all scaling factor). Figure [6] is a plot of these contours along with a colour map 



(2 + cos(27))e(l-0 = C/52, 



(51) 
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Figure 6. Contour plot of (I51I I illustrating the slow oscillations in amplitude and 
phase. The colour map shows the values of the contours. 



indicating the values of C/62. Notable features in this plot are the two maxima at 
(^,7) — (5,0) and (^,7r), where C/62 = 3/4, and the saddle point at (^,7) = (5, f ) 
where C/62 = 1/4. The maxima and saddle point are fixed points representing coupled 
motion where the amplitudes and phases do not change with time - stationary coupled 
modes. 

The two maxima correspond to motion along the diagonals of the guide. The 
oscillations in x and y have equal amplitudes and are either exactly in phase (7 = 0), 
or exactly in anti-phase (7 = tt). Even though the coupling is nonlinear, we have found 
the normal modes familiar from the theory of linearly coupled harmonic oscillators. 
Indeed, applying the multiple scales method to a pair of linearly coupled harmonic 
oscillators results in the relationship (1 + cos(27))^(l — — constant, which is 
remarkably similar to (|5ip . The frequencies of the normal modes can be found by 
differentiating their phase, 0(s), with respect to s. The phase of the X oscillation at 
position s is, according to (|45ll . <j){s) = ^s + 9i{cr), and so the oscillation frequency is 
oj = dip/ds — n + edOi/da. With the help of (|47cp the normal mode frequencies are 
found to be 



-NM =,[l- ^) . (52) 

We see that the frequency is lowered when x (and hence £62) is positive because the 
coupling term reduces the trap depth along the diagonals of the guide. The opposite 
is true when x is negative. Note that the normal modes with 7 = and 7 = tt 
have equal frequencies, unlike the linear coupling case where the coupling removes 
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the frequency degeneracy. There is another important difference too: the frequency 
of the normal mode differs from /i by an amount that depends on the square of the 
oscillation amplitude (the total energy). This is similar to the frequency shift caused 
by the cubic term discussed earlier, and quite different from a linear coupling whose 
normal mode frequencies are independent of amplitude. 

The two maxima in figure [6] are encircled by contours with 1/4 < < 3/4. 

In this range, the two transverse oscillations exchange energy completely. Looking 
at the contours near 7 = 0, the phase difference is zero when the amplitude of one 
mode has its maximum value, zero again when it has its minimum value, and reaches 
a maximum when the amplitudes of the two modes are equal. The phase difference 
never exceeds 7r/2. Again, this type of coupled motion is similar to that of linearly 
coupled oscillators. For a given value of 0/1)2 in this range, the maximum value of f 
is found at 7 = and is 



The saddle point in figure[B]is a coupled stationary mode corresponding to circular 
motion in the xy plane - we call it a circular mode. Its stationary nature is unstable in 
the sense that an infinitesimally small deviation from this point results in oscillations 
with time-varying amplitudes and phases. We will show later that the circular modes 
can become stable when 61 7^ 0. 

When C/62 < 1/4 the contours are below the saddle point and so no longer 
encircle the maxima. In this range, although the amplitude of one mode does grow at 
the expense of the other, the amplitudes do not exchange completely. The contours 
are open and the phase difference increases continually. Of particular interest is the 
fact that motion with ^ = or ^ = 1 is of this type, and does not yield any exchange 
of amplitude at all. The motions uncouple when all the energy is in one direction, 
in stark contrast to the classic motion of linearly coupled oscillators which exchange 
energy fully even when one oscillator begins at rest. For a given value of C/62 < 1/4, 
the large ^ contour has maxima at 7 = and tt, whose value, ^mi, is given by ([53|l . 
The small ^ contour has a maximum at 7 = 7r/2 with value 



Due to the exchange of energy between the transverse modes, there will be some 
particles whose maximal amplitudes of oscillation take them outside the boundaries 
of the guide, even though their initial amplitudes are within these boundaries. These 
particles are lost from the guide in the presence of coupling, and would not have 
been in its absence. We note again that the value of eh2 dictates the rate at which 
the amplitudes exchange, but not the degree of exchange. If the guide is infinitely 
long, and d)2 is small enough that our approximations hold, the loss of acceptance is 
independent of the size of e&2. A larger value only results in the particles being lost 
more quickly, but those same particles will eventually be lost even when the coupling 
term is infinitesimally small. Similarly, and somewhat counter-intuitively, the same 
particles are lost whether the nonlinear term increases or decreases the depth of the 
guide (along the diagonals). 

At the present level of approximation, it is straightforward to determine whether 
a particle with given initial conditions will be transmitted through the guide. We first 
relate the invariant quantity, C/62, to the initial conditions. To lowest order in e, the 
initial values of ai , 9i are related to the initial conditions via 




(53) 




(54) 
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ri(0)2 = X(0)2 + X'(0)V/i' 



(55) 




(56) 



and similarly for 02, ^2- With these, we can form ^(0) and 7(0) and hence C/62 from 
(|5ip. An equivalent approach is to realize that, to lowest order in e 



Knowing C/62 we use ((53)) and (|54p to find the maximum value of ^ for the particle, 
and hence the maximum values of ri and r2 ■ These maxima will eventually be reached 
(unless eb2 — exactly), and so the particle will be guided indefinitely if these are 
smaller than the boundaries of the guide, and will be lost with certainty if they are 
larger. 

We applied this algorithm to a large number of particles to find the set that are 
transmitted by a guide in the presence of a very small transverse coupling. We find that 
the coupling reduces the acceptance of an infinitely long guide to 72% of its acceptance 
in the absence of nonlinear forces. This is a general result that depends only on the 
coupling term being small enough that our lowest order approximation is valid. In the 
absence of couphng, the acceptance ellipses in the {X, X') and {Y, Y') phase spaces 
are uniformly dense - all particles that lie inside both ellipses are transmitted. This 
ceases to be true when the coupling is turned on. In fact, the sizes of the ellipses that 
enclose the stable particles do not change, but these particles are no longer distributed 
uniformly within them. For example, a particle that is close to the boundary of the 
ellipse in (X, X') may be driven out of the guide if it also happens to lie close to the 
boundary in (Y, Y'), but it will be successfully transmitted if it happens to lie close to 
the origin in (Y,Y'). Thus, the phase space acceptance in one direction has a dense 
region around the origin, becoming less dense towards the boundaries. 

By eliminating the independent variable, cr, from and we have obtained 
a great deal of insight into the coupled motion of the particles. However, we do not 
yet know the wavelength of the slow oscillations of amplitude and phase. To find this, 
we need to study the solutions of (|48l) and (|49)) with bi — 0. There appears to be no 
straightforward analytical solution, so we solved the equations numerically for a set of 
initial conditions corresponding to the full range of the contours seen in figure [6l and 
for several values of the single parameter appearing in these equations, e62i?o/(4/x) . 
In this way, we found the wavelength for the slow energy-exchanging oscillations (in 
units of s) to be 



Here, 5 is a number that depends only on the value of C/52, i.e. on the particular 
(^, 7) contour that the particle travels along. For most of the contours, the value of 
g is between 5 and 10. The shortest wavelengths {g ~ 5) are obtained for contours 
near the two maxima in figure [6l and near the edges of the plot where ^ is or 1. 
The longest wavelengths are found for values of C/b2 near 1/4 where the crossover 
from the closed contours to the open ones occurs. Taking an example where fj, = 0.5, 
eb2 — 0.1, Eq = r'j^ + r\ = 1 and g = 5 we find As to be 100 unit cells of the guide. 



C _ 3(^- 
b2 ^ 



.^XY + X'Y'f + fi^X'Y - XY')^ 



(57) 
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Figure 7. Variation of the phase space acceptance with the value of x- The guide 
has 100 lenses and has rjL = 1, rjS = 0. The points joined by a solid line show 
the results of numerical calculations with the micromotion included. The points 
joined by a dashed line were calculated in the small micromotion approximation. 



7.2. Larger coupling coefficient 

Having understood how the acceptance is reduced by a very small nonlinear coupling, 
we now consider larger values of the coefficient x- Applying the multiple scales method 
to higher orders in e becomes quite awkward mathematically, so we turn to numerical 
methods. Then we have a choice: we could either solve the exact differential equations 
given by (|24|) . or we could use the small- micromotion approximation and solve (|26p 
instead. In the former case, the step size needs to be chosen small relative to the 
micromotion wavelength, while for the latter it need only be small relative to the 
macromotion wavelength. For this reason, the approximate method is computationally 
faster by the ratio of the two wavelengths. Figure [7] shows the results obtained by 
both methods for a gapless guide containing 100 lenses each having rjL = 1. The two 
methods give similar results, the greatest difference being for large positive values of 
X where neglect of the micromotion underestimates the true acceptance. 

In line with our discussion in section [Y.li the acceptance drops on either side of 
the X — ^ position, due to the exchange of energy between transverse modes. This 
loss mechanism is fully established once x is large enough that a cycle of the slow 
amplitude and phase oscillations is completed within the length of the guide. The 
width of the sharp peak depends inversely on this length, a longer guide resulting in 
a narrower central peak. Away from the central peak there is a striking asymmetry 
between positive and negative values of x- This happens because there are two effects 
that cause the acceptance to depend on x, one which reduces the acceptance for both 
positive and negative Xj ^-nd another which reduces the acceptance on the positive 
side but increases it on the negative side. The first effect is the exchange of energy 
between the two transverse oscillations, already discussed above. Importantly, the 
wavelength of the amplitude and phase oscillations. As, is a function of the particle's 
transverse energy, Eq, and its value of C/62- Those particles with C/62 near 1/4 will 
have longer energy-exchanging wavelengths than others, and so this loss mechanism 
continues to be active as |x| increases. As a result the central peak in the figure has 
broad wings where there continues to be more particle loss with increasing |x| , reducing 
the acceptance on both sides of the central peak. The second effect is a direct result 
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Figure 8. Trajectory of a particle in the Jf Y-plane, for a gapless guide with 
riL = 1. The initial conditions were arbitrarily chosen. When x = the trajectory 
is a closed ellipse (thick line). When the transverse motions are coupled (thin line, 
X = —0.05), the trajectory is more complicated and the particle visits many more 
points in the plane. 

of the change in shape of the effective potential with the value of x- ^-S x increases, 
the nonlinear term makes the effective potential shallower along the diagonals of the 
guide, allowing an increasingly large fraction of the particles to escape along these 
directions. For negative values of x, the guide becomes deeper along the diagonals, 
allowing extra particles to be confined. It so happens that, for negative the two 
effects nearly cancel so that the acceptance is nearly a constant on the negative side. 

It is also important to note that the coupling term leaves the shape of the effective 
potential unaltered along the principal axes. For all negative the potential depth is 
shallowest along these axes. In the presence of the coupling, particles tend to explore 
many different regions of the guide, and so for many particles, increasing the potential 
depth along the diagonals does not improve the confinement since these particles can 
still escape along the principal axes. Figure [5] shows an example of the motion in the 
XK-plane when rjL = 1, comparing x = —0.05 to x = 0. In the absence of coupling 
the particle moves in a closed ellipse. The coupling results in a far more complicated 
trajectory, with the particle visiting far more points in the plane. 

Finally in this section, we consider the effect of the transverse coupling on the 
example guides that we examined earlier. In the case of the two-rod guide used in 
many previous experiments on focussing, guiding and decelerating, x is always negative 
and therefore has little effect on the acceptance. This, along with its simplicity, is a 
great advantage of the two-rod geometry. The four rod guide with two-pole symmetry 
discussed in reference jl^ has x = 0.19 when the Stark shift is linear, and this reduces 
the acceptance to about 40% of the value it takes in the absence of the coupling (for 
r]L = 1 and rjS = 0). When the Stark shift is quadratic x = 0.055 in this geometry 
and so the normalized acceptance is about 70%. The four rod bent guide used for 
guiding molecules in reference [11] has x = 1-5 when the Stark shift is linear. This 
very large value for x is a great disadvantage of this geometry since it reduces the 
acceptance enormously. It is, of course, necessary to consider the higher order terms 
to get a complete picture in this case. 
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8. Combined effect of both nonlinear terms 

8.1. Small nonlinear coefficients 

We return now to the general case where both bi and &2 are nonzero. Figure [S] shows 
contours of (|50|) for several different values of the ratio &1/62, along with colour maps 
to indicate the values, C/62, of the contours. As before, the fixed points in these plots 
correspond to the coupled stationary modes. The ones at (^,7) = (i,0),(i,7r) are 
the diagonal in-phase and anti-phase normal modes, while the fixed point at (5, f ) 
represents the clockwise and anticlockwise circular modes. When 61/62 is small, as in 
(a), the plot closely resembles the 61 = case shown in figure [51 The main difference 
is that a greater fraction of the contours encircle the two maxima that correspond to 
the normal modes. Part (b) of the figure shows the contours for 61/62 = 0.5. The 
saddle point has been replaced by a minimum meaning that the circular modes, as 
well as the diagonal modes, are stable. All motion involves a complete exchange of 
energy between one transverse direction and the other. As 61/62 is increased further, 
the area of the plot containing contours centred on grows. Eventually, once 

61/62 = 1 as in part (c) of the figure, this area includes all the contours. With a ratio 
greater than 1, e.g. 61/62 = 2 as shown in (d), the diagonal modes become unstable, 
and a larger fraction of the contours are associated with these modes. 

It is clear from figure [5] that the coupled stationary modes can be stable or 
unstable, depending on the ratio of 61 to 62. To find the stability condition, we 
form the Hessian matrix of the left-hand-side of ((50)) . A fixed point is a stable one if 
the determinant of this matrix is positive. Using this procedure it is straightforward 
to show that the normal modes are stable provided 61/62 < 1, and the circular modes 
are stable when 61/62 > 1/3. 

We can also easily find the spatial frequencies of the coupled stationary modes, 
following the same procedure as before. The frequency in the X direction is 
w = /i -I- edOi/da. With the help of (j47c|) we find the normal mode frequency to 
be 



Comparing these with equations p6)) and ((52)l . we see that, to this order in e, the 
change in frequency resulting from the inclusion of both cubic and coupling terms 
together is just the sum of the changes introduced by each term alone. 

8.2. Larger nonlinear coefficients 

To find the 2D phase space acceptance over a wide range of C, and x we solved the 
equations of motion for a large number of particles using the effective potential given 
by The results for a 100- lens gapless guide with r^L = 1 are shown by the 

points in figure [H for -0.025 < C < 0.15 and -0.3 < x < 0.3. For aU values of 
C, the results are rather similar to the C = result shown in figure [7] - there is a 
narrow spike centred at x = 0, a decline in acceptance for positive x and a plateau for 
negative x- This similarity suggests a hypothesis, which we could write symbolically 




(59) 



and the frequency of the circular mode to be 




(60) 
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Figure 9. Contour plot of H50| l illustrating the slow oscillations in amplitude and 
phase, (a) bx/b2 = 0.25, (b) fei/fe = 0.5, (c) fei/fe = 1, (d) 61/62 = 2. 

as A(C, x) = -^oP^iOPxix)- words, the hypothesis is that the total acceptance is 
the product of the one obtained from the hnear theory, Aq^ the factor Fq given by 
(j426p and (|42cp . which is independent of x and accounts for the cubic term, and a 
third factor, F^, that accounts for the couphng term independently of C, and is plotted 
in figure [T] The lines in figure [TO] are plots of this product form, showing that for a 
large region of the parameter space this simple method agrees with the full numerical 
calculation to within a few percent. The largest discrepancies occur when both x a-nd 
C are large and positive. Here, the numerical calculations show that the acceptance is 
considerably larger than the product form suggests, e.g. 70% larger when C, = 0.15, 
X = 0.3. Away from this region, the next largest discrepancies occur when x is small 
and negative. The numerical calculations show that in this region the acceptance 
tends to be larger than the product form predicts, by up to 20%. 

We are, at last, ready to put everything together to find the best way of building 
alternating gradient guides, decelerators and traps. We first found that, in the linear 
theory, the 2D phase-space acceptance is proportional to 77^, and so to the value of 
jasj. However, making 03 as large as possible may not be the best strategy, because of 
the detrimental effects of positive C and x which are related to one another through 03 : 
X + iC — 2?T.|a3|. If we make 03 too large, the increased acceptance suggested by the 
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Figure 10. Log-linear plot showing the 2D phase space acceptance of a guide with 
riL = 1, riS = 0, plotted as a function of x for various values of C,. The acceptance 
has been normalized to the value obtained when both x and C, are zero. The points 
are the results of numerical simulations using the effective potential, II29I I. The 
solid lines show the product of the numerical acceptance calculated with C, = 
(as in figure [Jll with the analytical result for the cubic term, equations ||426|| and 
l|42cll . 

linear theory will not be realizable due to the resulting large value of x- Decreasing x 
by increasing C, does not necessarily help, since large positive C, also reduces acceptance. 
In the effective potential picture, increasing 03 increases the depth of the potential, 
but increasing either C or x decreases the depth. We are only free to choose two of the 
three parameters, as, Q, x- What values should we choose to optimize the acceptance? 
We will look at some examples for particular values of the operating parameters, rjL 
and rjS. 

We start with the two-rod geometry. This is the easiest case to consider since there 
is only one free parameter, 03, determined by the the ratio of the rod radius, R, to the 
gap between the rods, 2ro: = [r^/ R)/{2 + tq/R) and 05 = a|. The two nonlinear 
coupling coefficients are therefore C — 2|a3| and x = — 4|a3| when the Stark shift is 
linear, and C = Sjasj, x — ~^\o-3\ when the Stark shift is quadratic. The calculation 
is hugely simplified by x being negative because we know that the acceptance is 
insensitive to x when it is negative, and we have an approximate analytical expression 
to describe the effect of the cubic nonlinearity. Our question can therefore be answered 
analytically with good accuracy. 

Figure fTTT a) shows the transverse acceptance versus 03 for a two-rod gapless guide 
with 77L = 1, in the case of a linear Stark shift. The unit of the acceptance is TZr^ 
and is independent of 03. See equation ([¥]) for the definition of TZ. The solid line in 
the figure is the analytical result obtained in the absence of nonlinearities, (|17p . The 
dashed line in the figure is also a completely analytical result - we take the result given 
by the solid line, multiply by the correction factor given by (|426p and (|42cP to account 
for the cubic nonlinearity, and finally multiply by 0.9 to account for the negative 
value of X, as suggested by figure [71 As 03 is increased from zero the acceptance 
first increases linearly, just as expected from the theory when the nonlinearities are 
neglected. However, as 03 increases further, the growing cubic nonlinearity halts and 
then reverses the ascent of the acceptance, which reaches its maximum value near 
03 = 0.03. We note that the rod radius needs to be 16ro to obtain this value of 
03. The maximum acceptance is approximately 0.0087?.rQ. We will examine this 
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Figure 11. 2D acceptance versus aa for an AG guide with a two-rod geometry, 
in the case of a hnear Stark shift, (a) Guide with no gaps: rjL = 1, rjS = 0. (b) 
Guide with gaps and lenses of equal length: rjL = r]S = 0.7. Dashed lines are 
analytical results. Round points are results of numerical simulations based on the 
effective potential. Square points are results of full numerical simulations. The 
solid line is the result obtained in the absence of nonlinearities. 



value in the context of some typical experiments at the end of the paper. The round 
points in the same figure show the results of numerical calculations of the transverse 
acceptance using the effective potential approach. Our analytical result agrees very 
well with these numerical calculations. The square points show the results obtained 
from a full numerical calculation, including the micromotion. As in previous examples, 
the inclusion of the micromotion increases the acceptance, but the overall shape of the 
plot remains unchanged, and the maximum of the acceptance still occurs close to the 
same value of 03. We have shown the result for a gapless guide operated at tjL = 1 
which is close to the optimal lens length (see figure [1]). Results for other cases are 
just as easily obtained using the same procedure. In general, the optimal value of 03 
increases as 7]L is increased towards the edge of stability. Figure fTlT b) provides the 
same information for a guide with gap lengths equal to the lens lengths. This is the 
situation typically found in AG stark decelerators. Again, we chose the parameters 
close to the optimal ones in the linear theory, specifically rjL = rjS = 0.7. The figure 
shows exactly the same trends as for the gapless case, but the acceptance is now 
maximized when 03 ~ 0.07 — 0.08, requiring rods of radii 6-7ro. Again, the numerical 
simulations based on the effective potential agree very well with the analytical result. 
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Figure 12. 2D acceptance as a function of as, wlien "qL = 1, riS = 0, calculated 
for several values of C,. The dashed lines simply join the points to improve the 
clarity of the figure. The solid line gives the result in the absence of nonlinearities. 
(a) Linear Stark shift, (b) Quadratic Stark shift. 



while the full numerical simulations show that the true acceptance is a little larger. 

Returning to the more general case, we will take as and C, to be the two free 
parameters, x being fixed once these two are chosen. We will concentrate on the 
example of a gapless guide operated with rjL = 1. Then, there are no new calculations 
to do since all the information we need is already given in figure 1101 For each value 
of aa and C, we calculate the acceptance in the linear case, which is just linear in 03, 
and then read off the multiplication factor from figure [TUl 

We present the results, for both linear and quadratic Stark shifts, in figure [121 
Let us concentrate first on the linear case, part (a) of the figure. We know from figure 
[4] that negative values of C, increase the acceptance above the value given by the linear 
theory, and it is tempting to conclude that we should design the guide with negative 
However, negative Q implies large positive values for x which greatly reduces the 
acceptance. As a result, the optimal acceptance is small and occurs for small values 
of 03, as illustrated by the results for C = —0.025 shown in figure flWa). The figure 
shows that increasing the value of C, increases the maximum acceptance, and also the 
optimal value of 03. This is because the beneficial effect of a smaller x outweighs 
the detrimental effect of a larger C. The acceptance reaches its maximum value of 
0.0257?,ro when C, ~ 0.075 and 03 ~ 0.11. At this point, x is very close to zero. 
For values of 03 below the optimum, x is negative and so has very little influence on 
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the acceptance, which consequently increases rather hnearly with a^. As soon as x 
becomes positive further increases in 03 only reduce the acceptance because of the 
associated increases in x- When C, > 0.075, the same trends are seen - the acceptance 
increases linearly with 03 up to the point where x is zero, and then decreases. Because 
increasing ( reduces the acceptance, the slope of this linear rise is smaller than before, 
the acceptance is maximized at higher values of 03, and this maximum value is reduced. 
For comparison, the result obtained in the absence of nonlinearities is shown by the 
solid line in the figure. 

The results for a quadratic Stark shift are shown in figure [TWb). The figure shows 
exactly the same characteristics as for the linear Stark shift, but with all features 
occurring at lower values of 03. The acceptance obtains its maximum value when 
C ~ 0.075 and 03 ~ 0.06. As before, x — at this point, and this is also true of the 
turning points for higher values of C,. The maximum acceptance is 0.0117?.rg, about 
half the value obtained for the linear Stark shift. 

The recipe for designing an AG guide should now be clear - find an electrode 
structure that produces a small positive value of ( (close to 0.075), and a x very close 
to zero, and choose lens lengths such that rjL is close to 1. If there is a need for gaps 
between the lenses, such as in an AG decelerator, it is best to keep the gaps short, in 
which case the above conclusions are not much altered. 

9. Conclusions 

In this paper we have shown how nonlinear forces alter the trajectories of neutral 
particles in an AG guide, and have shown how to calculate the phase space acceptance 
of guides with these nonlinear forces included. We have found how to design the guide 
so that the acceptance is optimized. To conclude, we apply our results to some example 
experiments, revisiting those discussed in Sec. 13.21 

Consider a beam of CaF molecules entering an AG guide with a forward speed 
of 350 m/s. The guide has ro = 1 mm, Eq = lOOkV/cm, r]L = 1 and 775" = 0, and the 
electrode geometry is designed so that 03 and C, have optimal values, as in figure fT2l 
The Stark shift will be linear in this case, and the ratio of the Stark shift at the centre 
of the guide to the forward kinetic energy is TZ = 1/88. The transverse acceptance 
will be approximately 0.0257?,rQ = (17mmmrad)^. The equivalent trap depth of this 
guide is approximately 30 mK. The available acceptance is well matched to the phase 
space area occupied by the beam in a typical supersonic beam experiment where the 
molecules pass through a small skimmer before entering the guide. We can therefore 
expect the guide to transmit a large fraction of the total beam. For a second example 
we consider Li atoms launched with a speed of 10 m/s from a MOT, into an AG guide. 
The guide has the same parameters as above, except that 03 is halved so that the 
acceptance is optimized for the quadratic Stark shift. In this case TZ = 0.47, and the 
transverse acceptance is approximately 0.017?.rQ — (70mmmrad)^. The equivalent 
trap depth of the guide is approximately 50 /iK, similar to the temperature of the 
atoms in the MOT. Again, we can expect the guide to transport a large fraction of the 
atoms. Unless the length of the guide is shorter than the macromotion wavelength, in 
which case the transmission will be larger than the above estimates, the acceptance 
does not depend on the length. The guide can be made arbitrarily long, delivering 
atoms and molecules to an experiment far removed from the source. Our examples 
demonstrate that a carefully designed AG guide has an acceptance that is adequate 
for many common experimental situations, making it a useful tool for cold atom and 
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molecule physics. 
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